Correlation functions of the BC Calogero— Sutherland model 



Shinsuke M. Nishigaki,^'^ Dimitri M. Gangardt,^'0 and Alex Kamenev^'l^ 

'Department of Physics, University of Connecticut, Storrs, CT 06269, USA 
^LKB, Laboratoire de Physique de I'ENS, 24 rue Lhomond, 75005 Pans, France 
'^Department of Physics, University of Minnesota, Minneapolis, MN 55455, USA 

(Dated: September 4, 2002) 

fSJ ' The BC-type Calogero-Sutherland model (CSM) is an integrable extension of the ordinary A- 

' type CSM that possesses a reflection symmetry point. The BC-CSM is related to the chiral classes 

, of random matrix ensembles (RMEs) in exactly the same way as the A-CSM is related to the 

' Dyson classes. We first develop the fermionic replica a-model formalism suitable to treat all chiral 

^ ^ RMEs. By exploiting "generalized color-flavor transformation" we then extend the method to flnd 

^ , the exact asymptotics of the BC-CSM density proflle. Consistency of our result with the c = 1 

■ Gaussian conformal field theory description is verified. The emerging Friedel oscillations structure 

and sum rules are discussed in details. We also compute the distribution of the particle nearest to 
the reflection point. 
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I. INTRODUCTION 



It is a well known fact that there is an intimate relation between the one-dimensional quantum problem with 
the inverse-square interaction potential, i.e. Calogero-Sutherland model (CSM) j^], and Dyson random matrix 
CO . ensembles (RMEs) On the most elementary level the correspondence goes as follows: for the three particular values 
K>" ' of the coupling constant (A = 1/2, 1,2) square of the many-body ground-state wave function of the CSM coincides 
with the joint probability distribution (JPD) of RMEs with Dyson index fi = 2A. Consequently the knowledge of 
the RME correlation functions may be immediately translated to the information about the CSM. Historically this 
correspondence proved to be very fruitful for advancing the understanding of the CSM. 



. It was later realized that Dyson's classification of the RMEs was not exhaustive. Studies of the two-sublattice model. 



ps) I mesoscopic transport, and QCD Dirac spectra initiated the introduction of their "chiral" counterparts 
Subsequently, based on Cartan's classification of Riemannian symmetric spaces, Altland and Zirnbauer 
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have 



further added six "superconducting" chiral symmetry classes. These chiral classes are characterized by the special 
role played by the zero energy. Namely, there is a certain number of eigenvalues that must have zero energy, while 
all other eigenvalues occur in symmetric pairs (mirror images) around zero. As a result, the mean density of states 
(DoS) exhibits either a hollow or a bump around the zero energy followed by decaying oscillations at larger distances. 
Such a structure on the level of the mean DoS does not show up in the Dyson classes and, as we shall explain below, 
Q I may be called Friedel oscillations. 
^ , The question is whether one can find an appropriate generalization of the CSM whose ground-state wave function 
possess the same reflection symmetry property as the JPD of chiral ensembles. The answer is known to be affirmative. 
Indeed, one can write down an integrable one-dimensional model with inverse-square interaction, reflection symmetry 
^ ■ and special single-particle potential centered at zero having the required ground state. It is known as the BC-type 
' CSM in the literature ||l2|. Due to the presence of the mirror boundary and localized single-particle potential (an 
- - < impurity) the model lacks translational invariance and the resulting ground-state density is not uniform. In particular 
the density profile develops the Friedel oscillations far enough from the impurity. Accordingly one has a unique example 
of integrable strongly interacting models that exhibit Friedel phenomena. That gives one a possibility to gain the 
exact information on the amplitude decay rate, spectral characteristics and phase shifts of the Friedel oscillations in 
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the interacting system. Despite the proven integrabihty, exphcit form of the correlation functions of the BC-CSM 
with generic coupUng constant were not estabhshed, except for some partial results In], |l^, |l^ . 

The purpose of this paper it to fill this gap for rational values of the coupling constant. To this end we employ the 
recently developed approach based on the replica trick |l^, It was previously tested on the pair correlation 
function of the ordinary A-CSM where it perfectly agrees with the exact results of Haldane |20 and Ha ||2l| for 
any rational coupling constant. The idea is to explore the relation with the RMEs, where the replica trick was found 
to be accurate in the asymptotic regime. We thus develop first the fermionic replica approach to the chiral symmetry 
classes of RMEs. Not surprisingly, we are able to reproduce the asymptotic behavior of the known DoS profiles for 
chiral RMEs. We then extend the treatment away from the RMEs values of the coupling constant and obtain closed 
analytic results for any rational two-body coupling constant A and any impurity phase shift. 

We found that for rational values of A = p/q {p, q coprime), the spectrum of the Friedel oscillations contains exactly 
p harmonics, corresponding to 2kp^Akp, . . . ,2pkp density oscillations. The Z-th harmonic {I = l,...,p) decays 
algebraically as 9~ ' , where 9 is the distance from the impurity. The amplitudes of the harmonics depend on the 
number, Z, and the coupling constant, A, but are not sensitive to the strength and details of the impurity expressed 
through a phase shift, v. Moreover, the harmonics amplitudes are closely related to those of the two-point correlation 
function of the homogenous A-CSM. We provide a conformal field theory account for this fact. The impurity phase 
shift, V, affects only phases of the harmonics in the asymptotic regime. We explicitly show that the relation between 
the total charge expelled by the impurity and the asymptotic phase shifts, known as the Friedel sum rule, holds for 
the interacting system. We also compute the distribution of the locus of the particle nearest to the mirror boundary 
point for the BC-CSM at generic values of the coupling constants. 

The paper is organized as follows: in section || we briefiy introduce chiral ensembles, and develop the appropriate 
fermionic replica cr-models. Section ID is devoted to the introduction and replica treatment of the BC-CSM. In 
section IV we perform the analytical continuation and the replica limit, and extract the density profile. In section ^ 
we compare this result with the effective conformal field theory descripti on. In section ^ we compute the nearest 
particle distribution. A summary and discussions are provided in section VII. Technicalities of the calculations are 
relegated to the two appendices. 



II. CHIRAL CIRCULAR ENSEMBLES 



A circular RME is defined as an ensemble of unitary matrices U representing a Riemannian symmetric space P, 
stochastically distributed according to the Haar measure dU of V. Expressing a symmetric space as a coset V = G/H 
with a compact Lie group G and H <Z G, the Cartan mapping G — > G/H, g i-^- U{g) is as shown in the table below: 



class 


G 


H 


U{9) 


A (CUE) 


ViN) 




1 


9 


AI (COE) 


IJ{N) 




0{N) 


T 

9 9 


All (CSE) 


V{2N) 




Sp{2N) 


9^J9 


AIII (chCUE) 


V{N + N') 




V{N) X U(7V') 


I9^l9 


BDI (chCOE) 


SO{N + N') 




SO(iV) X SO{N') 


I9^l9 


CII (chCSE) 


Sp{2{N + N')) 




Sp{2N) X Sp(2iV') 


ig^ig 


D,B 


SO(2Ar),SO(2iVH 




1 


9 


C 


Sp(2iV) 




1 


9 


CI 


Sp(4iV) 




\Ji2N) 


ig^ig 


DIIIe,o 


SO(4iV),SO(4iVH 




U(2A^),U(2iV + l) 


9''9 



Here 



J = 


~1n 


, 1 = 


In 




_tN 




-1n' 



(2) 



{N 2N, N' 2N' for CII and A^, N' 2N for CI), and = Jg'^J'^ denotes the quaternion dual of the matrix g. 
Out of these twelve classes, last non-classical nine possess chirality |22j, i.e. non-zero cigcnphases appear in complex 
conjugate pairs. The JPD of these non-zero eigenphases 



U = l^diag(e'^i , . . . , e**" , e''^^ , . . . , e"**" , 1, . . . , l)y^ 



(3) 



for the circular ensemble is given by psl {0 < < tt) 



N 



PiOi, ...,9N)dei...deN = n sin^' T T I^Jv(cos 



1=1 

N 



(4) 
(5) 



i=l 



where A^r is the Vandermonde determinant of the rank N and the constants ci , C2 and (3 for all nine chiral ensembles 
are given in the table below. 



class 




Cl 


C2 


P 


V 


V (circular RME) 


M (F-rephca NLcrM) 


AIII 


2{N' - 


-iV)4 


- 1 1 


2 


N' - N 


V{N + N')/ (JJ{N) X U(A^')) 


U(n) 


BDI 


N 


-N 





1 


N' - N 


SO(iV + TV')/ (SO(iV) X SO(iV')) 


U(2n)/Sp(2n) 


CII 


4(iV' - 


-iV)4 


-3 3 


4 


N' - N 


Sp(2(7V + N'))/ (Sp(2iV) X Sp(2iV')) 


U(2n)/0(2n) 


D 










2 


-1/2 


S0(2iV) 


S0(2n)/U(n) 


B 




2 





2 


1/2 


S0(2iV+ 1) 


S0(2n)/U(n) 


C 




2 


2 


2 


1/2 


Sp(2iV) 


Sp(2n)/U(n) 


CI 




1 


1 


1 


1 


Sp(27V)/U(7V) 


Sp(2n) 


Dllle 




1 


1 


4 


-1/2 


SO(4iV)/U(2Ar) 


S0(2n) 


DIII„ 




5 


1 


4 


1/2 


SO(4iV + 2)/U(27V+ 1) 


S0(2n) 



(6) 



The JPD Eq. (j^) defined over y G [0, 1] may be comprehensively called Jacobi ensemble. The constant v, defined as 

Cl + 1 



f3 



-1, 



(7) 



may be called "topological charge" , when the tangent element of a random matrix U in the first three cases is 
interpreted as the QCD Dirac operator in even dimensions 

The real characteristic polynomial, or so called the fermionic replicated partition function, for these circular en- 
sembles is defined as 

Zn,N{0) ^ j dU Aei{S - e~'^U)" = j dg det(e*^ - e-'iu{g)Y . (8) 

V G 

After the color-flavor transformation and the thermodynamic limit where iV — > oo, ^ with their product 
(denoted by the same 9 for the sake of simplicity) fixed finite, it takes the form. 



Zn{9) = lim Zn,N , 



du{A'sXuY%'^ 



(9) 



where M. are the 'dual' symmetric spaces of the fermionic nonlinear cr-models which are listed in the table (|g). 
Hereafter we suppress irrelevant normalization constants, that goes to unity in the replica limit, n — s- 0. The derivation 
of Eq. from Eq. (H) for the AIII, BDI, and CII classes is summarized in Appendix In order to derive Eq. (||) 
for symmetry classes whose pertinent color-flavor transformations are not immediately available in the literature, 
one could as well adopt an alternative method of magnifying the origin of the circular ensembles first (i.e. employing 
Gaussian ensembles), performing Hubbard-Stratonovich transformation and taking the thermodynamic limit. For a 
Gaussian ensemble treatment of the BDI and CII classes, see Ref. [ p^ . 

Performing the integration over angular degrees of freedom v oiu = wdiag(e*'^")tit, one obtains 



Zn(9) = 0" 



n 



iyq>a.-\-iQ cos t; 



n 



sm 



4//3 



(10) 



This expression is valid for all nine chiral symmetry classes. Note that Eq. (10) depends only on /3 and c\ but not 
on C2, because we have magnified the vicinity of the origin, = 0. Accordingly the class C gives the same Zn{9) as 
B, and the class CI the same as BDI at ci = 1, reducing nine symmetry classes of chiral RMEs to seven universality 
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classes. The number could be further reduced by introducing three universahty classes of Laguerre ensembles, having 
continuous v and discrete /? = 2, 1, 4. 

Consequently we succeeded in expressing the integration over the initial iV-variable JPD, Eq. (Q), through the n— 
fold integral, Eq. (10). We shall discuss its evaluation, analytical continuation and the replica limit after we introduce 
the BC-CSM. This replica treatment was previously performed for the AIII class in Ref. pa]. 



III. BC CALOGERO-SUTHERLAND MODEL 



The generalized Calogero-Suthcrland Hamiltonian |12 is defined as 



H 



V — 



1=1 



^ sin2(6l-a/2) 



(11) 



a6A 



Here A is a root system of a Lie algebra in A'^-dimensional vector space, d is the vector (^i, . . . , ^at), and g\a\ is a 
coupling constant depending only on the root length. Quantum integrability is ensured by these conditions. The 
ordinary, translationally invariant model corresponds to the Aat-i root system. 

The quantum one-dimensional model of N interacting particles on a semicircle, < 0i < tt, where i = 1, ... ,7V 
with the Hamiltonian 



H 



N Q2 



1(1^1) 



4=1 



N 

E 



1 



N 



^l(^l-l) f^ 1 



^2(^2-1) 



1 



. , sm ^ 



N 

^ cos^ ^ 

t—1 2 



(12) 



corresponds to the BCn root system, hence called the BC-CSM. As the BCat root system contains roots of length 
l,\/2,2, there are three independent coupling constants. A, Ai,A2. This family contains CSMs corresponding to Bjv 
(A2 ~ 0), Cat (Ai — A2), and Djv (Ai = A2 = 0) root systems as its subfamilies. In addition to the pairwise inverse- 
square interactions, the particles interact with their own mirror images (occupying the other semicircle, t: < 9 < 27r) 
and with the two single-particle impurity potentials placed at ^ = and 9 ~ n. In what follows we shall assume the 
thermodynamic limit, N — )■ 00, and shall focus on the vicinity of the impurity ai 9 = 0. The other impurity may be 
treated in exactly the same manner. The model is known to have the following ground-state energy 



N 



H^o^Eo^o, £^0 = ^ ((^-«)A + 

i— 

and the ground-state wave function 



Ai + A2 



N 



(13) 

(14) 
through 



The absolute square of the ground-state wave function coincides with the JPD of the chiral RMEs (^ 
a change of the coupling constants 

I ^[3/2, h=ci/2, ^2-02/2. (15) 

Notice, that Eq. (|lj) is not restricted to the special values of /3, ci,C2 listed in the table (^. The energy and wave 
functions of the excited states were studied in Refs. |30[ |. 
The particle density is defined as. 



N 



N 



{p{9)) ^{Y,5{9-9,)) = v/^(T^ ( E ^(y - y^) 



(16) 



where y = siv?{9/2). The angular brackets denote ground-state expectation values, or equivalently averaging over 
the normalized JPD, Eqs. (^ and (|^), for the first and second equalities correspondingly. One may then employ the 
replica trick to write 



N 



1 d ^ 
lim — 5to— JJ(y — Vj — ie)" ■ 



n^O riTT 



dy 



(17) 
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As a result, one obtains 



n-*o rnr ay 



(18) 



j,=sin2(6l/2) 



where the "repHcated partition function" is defined as 



r " 

ZnMy) = y n {dy^ vt'^'i^ - v^f'-"\y ~ yO") {^N{yfy 



(19) 



Baker and Forrester |Q have noticed the integral equality due to Kaneko and Yan , which we suggestively 
call the "generalized color-flavor transformation" . With its help one may express the partition function in the following 
way 



f " f 'l+'2 + l o 

Zn,N{y)= [[\dXaXa ' {I - Xa) 



l2+n-l/2 



Xa{l - yxg) 
I - Xa 



(20) 



The integration contour C encircles the cut between Xa = and = 1- The general form of the integral identity is 
given in Appendix So far no approximation has been made. Now we pass to the thermodynamic limit, — s- oo, 
and magnify the vicinity of the 9 — impurity. To this end we rescale the variable as — > 9/N and correspondingly 
y ~ 9"^ /{AN"^). By redefining the integration variables as Xa = 1 — 2iA6'~^e*'^" and taking the thermodynamic limit, 
N —I- oo, one finds (see Appendix IbI for details) 



Zn{6) = hm Zn,N 



N 



)« n 
i a=l 



n 

.a>b 



sm 



l/l 



One may introduce notation 



A 2A 



(21) 



(22) 



to notice the exact coincidence with the a-model representation of the chiral RMEs, Eq. (10), provided that the 
coupling constants are related via Eq. ([l^) . The important difference is that the BC-CSM representation in the form 
of Eq. (pl|) is not restricted to the RMEs values A — 1/2, 1, 2 and special values of the topological charge, v. 



IV. ANALYTIC CONTINUATION AND REPLICA LIMIT 



Consider the n-fold integral 



y n 

Zn{9)^d''- J Hid^a' 



iu4'a,~\'iO COS 0a 



n 

.a>b 



sm 



<t>b 



l/A 



(23) 



One may stretch the integration contour from the unit circle in the complex plane of Za = e*"^" into two lines parallel 
to the imaginary axis with 3ffe — ±1. The original integral, Eq. (p3|), splits into the sum of n terms with I integrals 
having 3?e Zq = — 1 and remaining n — I ones 3?e = 1; here I — 1, . . . , n. The further progress is made possible in the 
asymptotic limit, 3> 1. In this case the integrals are dominated by the vicinities of the saddle points z^ = ±1 and 
thus may be evaluated employing the Selberg integral. This strategy was described in details in Refs. |jl^, |l^, p^ . 
Proceeding this way, one finds for the replicated partition function 



'^F^(A) e"'^'-2»ie2- 



(24) 



where we have omitted a normalization constant that goes to unity in the replica limit, n 
1^, |l^, |l|, |l^ , introduced the notation 



and, following Refs. 



(25) 
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Employing the observation that -F,j(A) = for / > n, one may extend summation over I in Eq, 
then perform the analytic continuation n — + 0. As a result, one finds via {p{9)) — lim„_,o(7rn 



_J) to infinity and 
^'imdgZnie) (cf. Eq. 



1 



1 + 2 



2 2A 



(26) 



where 



2iyx lir(l-(a-l)/A) 



(27) 



Eq. ( |2q ) for the asymptotic of the ground-state density of the BC-CSM is the central result of this paper. For the 
integer values of the coupling constant, A, it may be extracted from the expressions derived by Baker and Forrester 
|l|, see also Ref. |l|. 

We remark that each term in Eq. (|2^ ) is the one with the lowest power in 6^^ among all terms carrying the same 
frequency 21. This corresponds to truncating all contributions from the descendent fields in the conformal description, 
see the next section. One could compute these secondary terms, subleading in 9, by performing a perturbative 
expansion around each saddle point with the help of the loop equations Eq]. 



V. CONFORMAL FIELD THEORY DESCRIPTION 



By comparing the one-particle correlation function of the BC-CSM (^6|) with the equal-time two-particle correlation 
function of the ordinary A-CSM M, |2l 



(P(^)P(O)) 



1 



{2n) 



1 



2V^^cos(2;0) 



2A6'2 ' e^^y^ 
1=1 



(28) 



one notices that it is expressed through the very same coefficients dz(A), Eq. (|27|). It is thus clear that the harmonics 
amplitudes, di{X), are properties of the homogenous interacting system and not of the localized impurity. As explained 
below, this result could be anticipated from the effective conformal field theory description. The latter is capable of 
predicting the low energy properties of the system apart from numerical values of the coefficients. 

Based on the finite-size scaling analysis, Kawakami and Yang |34 identified the low-energy effective theory of 
the CSM in the thermodynamic limit to be the c = 1 Gaussian conformal field theory at radius R — VA/2, either 
non-chiral (A-CSM) 



zvr 



$(z,z) 



(z) , $ = $ + 2txR , 



(29) 



or chiral (BC-CSM) |lj, Namely, they have found the density operator should have zero winding number. It 
does not have a definite conformal dimension, and therefore is expanded in terms of primary and secondary operators 
whose left- and right-moving vertex momenta are equal. Here we shall consider only contributions from the U(l) 
current and the primary fields (vertex operators with charges allowed by the compactification (p9|)). 



p(z,z) 



Po 



6(a,0(z) + 9j0(z)) 



di e 



il{z+z) a(t){z)/R a(j>{z)/R 



(30) 



where the expansion coefficients, h and di[— d^i), are not determined from the conformal field theory. Neither are the 
oscillation factors e»'(^+^\ which are set by hand to describe transport of I pseudo-particles ||2l[ from the left Fermi 
point {—kp = —1 by normalization) to the right Fermi point (kp — 1) |3^. We have factored out po such that the 
constant term carries do = I. The propagator and the vertex correlator are given by 



/ (z - z'yyi^"'') 



(31) 



The coefficient in the second equation is a matter of convention and reflects a particular choice of the ultraviolet 
regularization. Another choice of the regularization would change coefficients di, but not the final result. Employing 
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Eqs. (30), (|3l|), one obtains for the equal-time two-particle correlation function of the A-CSM (we denote z 



-it) 



OO 

E 

^— — OO 



^2il(9-e') 



(32) 



Comparing this expression with Eq. (28), one finds that the b = l/\/A |2^, while di = di{X), cf. Eq. (27). Let us 
consider now the BC-CSM and concentrate on the case without the phase shift for simplicity. The Dirichlet boundary 
condition at die z = is translated into the open boundary bosonization rule p5| 



(z) for 5to z < 0. 



(33) 



As the right mover is identified as (—1) times the left mover at the mirror-imaged point, there exists a nonzero matrix 
element between the left- and right-moving vertex operators: 



Sw 



I [z + z'Yli^«') 
As a result, the mean density of the BC-CSM becomes nontrivial, 



{p{e + »o)) ^p^Y, di 



(26»)'V(4fl=) 



(34) 



(35) 



It is in exact agreement with Eq. i^x^ if one disregards the phase shift. The phase shift may also be included in 
the conformal description by shifting the identification of the right and left movers, Eq. (p3|), by a constant factor: 
-7r(i^ + l/2- 1/(2A)). 

We found that the knowledge of the asymptotic behavior of the correlator of the homogenous A-CSM, Eq. ([2^), 
supplemented by the conformal field theory description is, in principle, sufficient to predict the BC-CSM correlation 
function, Eq. (Pq). This agreement indicates the fact that the ultraviolet property of the field that is responsible for 
the normalization of vertex operators are not affected by the presence or absence of the boundary. Well anticipated as 
it is, we nevertheless consider this fact to be worth verifying, as done in this paper. This fact can be put on a further 
test by computing the asymptotics of e.g. two-particle correlation function for the BC-CSM, though it is technically 
more challenging. We could as well reverse the logic and conjecture the asymptotically expanded form of any p-point 
correlation function of density operators for the A-CSM {p > 3) or for the BC-CSM (p > 2), by using Eqs. (|30|), (IsT 
and M). 



VI. NEAREST PARTICLE DISTRIBUTION 



Another quantity of interest in the theory of interacting electrons is the probability E[s' , s] of finding no particle 
within an interval [s', s] or the particle spacing distribution p[s) that is a derivative of the former. In the context 
of spin chains a similar quantity was recently discussed in Ref. |]36|| . For nonchiral as well as chiral RMEs, Tracy 
and Widom |37j has developed Mehta's computation of E[s' , s] as a Fredholm determinant into a systematic and 
powerful method. As their method determines E[s',s] as a solution (r function) to a transcendental equation of 
Painleve type relies upon the orthogonal polynomials, its validity is necessarily limited to A = 1/2,1,2. On the 
other hand, E{s) = E[Q,s] for the chiral RMEs has been computed by an alternative and far simpler "shifting" 
method |^, ^ as explained below. We show that with a help of the generalized color-flavor transformation, 

this method is applicable also to the BC-CSM at generic values of the coupling constants. Consider for a moment 
Ai = n + 1/2, where n = 0, 1, 2, ... is an integer and A2 = 1/2. As we are interested in the universal behavior in the 
vicinity of the reflection point 6 = 0, the restriction on A2 is irrelevant. The probability of having no particle within 
an interval 0<6'<s, orO<y<F with Y = sin^(s/2), is defined as 

1 TV 

En{s)= const. ff[{dy,y^) {AN{y)^)^ • (36) 

The constant should be chosen to ensure Epf{0) = 1. By shifting and rescaling the integration variable as y ^ 
(1 — Y)y + Y, one obtains 



En{s) = const.(l - / J| 



Y ^ "1 



{^N{y)T ■ (37) 
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Now we apply the generalized color-flavor transformation (gO|) to reexpress E^^s) as an n-fold integral, 

N' 



/II, 



dXa Xa (1 - Xa) " 



1 



(A„(x)2)i/\ (38) 



We finally rescale s s/N ; Y s^/(4A^^) and take the thermodynamic limit. Following the same procedure that 
lead from Eq. to Eq. (|l|), one finds 



E{s) = lim En (^) 



const, e 4«^s"(i 



i+S COS (, 



.a>b 



4>b 



111 



(39) 



This exact result has previously been derived from Laguerre (chiral Gaussian) ensembles at arbirtary A and at 
Ai =n + l/2 @. 

To cornpute the asymptotics of E(^s) for s ^ 1, one can relax the restriction on Ai, by first evaluating the n~fold 
integral (|39| ) by the saddle point method and then performing the analytical continuation n — > Ai — 1/2. In the 
large-s limit we pick only the contribution of the replica-symmetric saddle point 0^ = 0, since the contribution of all 
other saddle points is exponentially smaller (the difference with the previous computation is that s enters Eq. ( |39| ) 
without imaginary unit). This way Forrester |0 has derived the asymptotic (s ^ 1) result 



E{s) 



2A e 4 e 



const, s 2 
const. 5-K^l-i)c-T^■'+(^l-4)^ 



(40) 



where v is the topological charge defined by Eq. (|2^). The Gaussian factor exp(— As^/4) could be anticipated from 
the mean-field treatment of the classical logarithmic gas . The other factors in the asymptotic expression Eq. ( |40| ) 
could not be found in any simpler way, to the best of our knowledge. 

The distribution p(s) of the locus s of the particle nearest to the reflection point is given by p(s) = —dsE{s). In 
the other limiting case, s ^ 1, p{s) is determined by the interaction of the particle closest to the reflection point to 
its own mirror image. Inspecting Eq. (^ one immediately finds p{s) cx s^^^ 



VII. DISCUSSIONS 

Let us now take a closer look at our main result, Eq. (p6|). The constant term on its r.h.s., po ~ I/tt, represents the 
uniform density of particles {N particles within [0, tt]) far away from the impurity. It may be traced back to the replica 
symmetric contribution (all n integrals are taken at = — 1 saddle point) to the partition function. The decaying 
(as 6 — > oo) oscillatory terms on the r.h.s. of Eq. (|2^) are the Friedel oscillations of the particle density induced 
by the mirror boundary and the impurity potential. These terms may be identified as the replica symmetry broken 
contributions to the partition functions (/ integrals are taken at the "wrong" saddle point, Za = 1). In general, there 
is an infinite number of harmonics (unlike a single "2fcj7" harmonic in the non-interacting system!) in the oscillation 
spectrum. Since di cx exp(A^^/^ Inl) for I ^ 1, the sum over harmonics on the r.h.s. of Eq. ( p6| ) is, in general, not 
convergent. It is not clear to us at the moment, whether there is a consistent resummation scheme. There is, however, 
an important class of the parameters, where Eq. ( p6| ) is mathematically rigorous. 

For any rational coupling constant A = p/q the coefficient di{p/q) = for I > p and therefore the sum terminates 
after exactly p oscillatory components. One finds that, in addition to the usual 2kF Friedel oscillation, the system 
possess Akp, . . . , 2pkp oscillatory components of the density (in the unit kp — 1 accepted here). This fact might be 
expected from the form of the density-density response function of the homogenous A-CSM However, the 

algebraic decay rate of the harmonics could not be determined employing linear response of the A-CSM. Indeed, the 
latter predicts that the l-th harmonic decays as 5/-2'~/a+i^ while the correct decay rate is 6*^' Notice, that for 
non-interacting particles, A = 1 and thus I — 1, both ways give the correct one-dimensional decay of the 2fci? Friedel 
oscillations: 9^^. For any interacting system, A 7^ 1, the linear response is bound to fail in the asymptotic regime. 
These observations was already made in the Tomonaga-Luttinger liquid literature ji^, ^ . Now we can confirm them 
having the exactly solvable model system. 

In the BC-CSM one finds not only the decay law, but also the relative amplitudes of the harmonics: coefficients 
di{l), Eq. (p7|). Notice that these amplitudes are determined by the interaction strength. A, only and are independent 
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on the impurity strength, Ai. This is due to the fact that the mirror boundary condition induces oscillations of 
the maximal possible amplitude. The additional single-particle potential centered ai 9 — and characterized by Ai 
changes the phase of the oscillations only. The entire information about the impurity strength, Ai, is incorporated in 



the parameter Eq. (22). In the asymptotic regime the latter affects the phase of the Friedel oscillations only and 
therefore may be associated with the impurity phase shift. (Unlike the leading order, the amplitudes of sub-leading 
perturbative corrections in negative powers of 6 do depend on the phase shift, v.) 

To verify Eqs. (p6|), ( p7| ) one may compare them with the available exact DoS of the chiral RMEs (see Ref. 
and references therein). Employing Hankel's asymptotic expansion of the Bessel function, one may check that these 
asymptotic perfectly agree with Eqs. (p^), (|27|). One may also notice that for D, B, and C symmetry classes we 
have obtained the exact rather than the asymptotic results. This coincidence is due to the Duistermaat-Heckman 
localization theorem Q. We see that having unitary symmetry class, /3 = 2, is not sufficient to satisfy the 
Duistermaat-Heckman theorem. One should also have the special value of the topological charge, v — ±1/2, to 
secure cancellation of all higher-order perturbative corrections . 

One may define the total charge attracted (expelled) by the impurity to (from) the region near ^? = as: 

oo 

Q = I de {p{e) - po) , (41) 



where po — 1/t: is the uniform asymptotic density. For the chiral RMEs, where the exact expressions including small 
9 region are available, the result is (cf. Eq. ( p6| ) ): 

^2V 2 2Ay4 2A ^ ' 

Notice, that the attracted charge depends both on the impurity amplitude, Ai and the interaction strength. A, reflecting 
the fact that the impurity is screened due to the interactions. The pure mirror boundary, without the single-particle 
potential attracts quarter of a particle irrespective to the interaction strength. Equation ( ^ ) is a manifestation of 
the famous Friedel sum rule: the total expelled charge is equal to the impurity phase shift (divided by 2tt); the latter 
also determines the phase of the density oscillations far from the impurity. We conjecture, in accordance with the 
earlier works Esl, that Eq. (E3) is valid for any values of I and Ai. 
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APPENDIX A: a-MODEL DERIVATION VIA COLOR-FLAVOR TRANSFORMATION 



We first consider the simplest case AIII with v ~ Q. The Grassmannian \]{2N)/ {\]{N) x \]{N)) is a complex 
Kahler manifold, and its unitary matrix representative U in Table (^ is conveniently parameterized by the complex 
stereographic coordinate Z.y , i, j = 1, . . . , A'^, as 



1 Z 

-Zt 1 



z e C 



NxN 



U = Ig^lg — ^, 7 = 

The Kahler potential k{Z, Z) ~ tr log(l + ZZ^) leads to the Haar measure 

dU = 

The replicated partition fimction then reads 



/ = 



1 

-1 



det(l + ZZt) 



2N 



(AI) 



(A2) 



Zn,N{9) — 



V{2N)/(V(N)y.V(N)) 



dU det(e'5 -e-'^uy^ 
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Ud'Z 



civxiv det(l + ZZt)2JV V I I J 

f H'^'^ 



(A3) 

We introduce a set of (A^xn)-coniponent independent Grassmannian numbers -0", xl^ i^ii xfi i — ^, ■ ■ ■ , N,a — 1, . . . ,n 
to exponentiate the determinant, {I'jl ~ 7^) 



(A4) 



Now we employ Zirnbauer's color-flavor transformation |2S 

Ud'Z 



C'NxN det(l + ZZt) 



2N+n 



exp I 



UtZ^.X^ - Xizy;) = / du exp (V;>°V^' + X^^^'X?) , (A5) 



to obtain 



dil;dtpdxdx / du exp (^(e'2 - e~'^)('il;ip + xx) + (e'^ + e~'^){ipuilj + xw^x)j 



U(n) 

N 



du det ( cos 61 + 1 sin 9 ) . (A6) 



U(n) 

We stress that no approximation has been made in the above procedure. 

In the thermodynamic limit, N 00 and — > 0, one has cos0 ~ 1 and sinS ~ 0. As a result, the determinant in 
Eq. ( [A^ ) may be exponentiated: 

Znie) = lim Zn,N ( ^) = [ due*^'^ (A7) 

U(n) 



N^oc ' \ N 



One can repeat the above procedure for the BDI and CII classes. The parametrization (Al) of the real and 
quaternionic Grassmannian manifolds, S0(2iV)/(S0(iV) x SO{N)) and Sp(4iV)/(Sp(2iV) x Sp(2iV)), involves x iV 
real and quaternion-real matrices Z , instead of complex. The color-flavor transformation trades the integrations over 
these 'colored' variables Z with the ones over 'flavored' variables u that are antisymmetric and symmetric unitary 
matrices, respectively (straightforward as they are, such types of color-flavor transformation have yet to be exhibited 
explicitly in the literature, to the best of our knowledge). Accordingly, the integration domain of the transformed 
partition function ( |A6| ) becomes antisymmetric unitary matrices (U(2n)/Sp(2n)) and symmetric unitary matrices 
(U(2n)/0(2n)), with the rest being unaltered. Inclusion of non-zero v is straightforward by consi derin g a rectangular 
Z. It merely shifts the the power of the determinant in the measure ( [A^ ) by ly, and modifies Eq. (A7) into Eq. (||) in 
the thermodynamic limit. 

For the B-D and C or Dili and CI classes, one utilizes the color-flavor transformation between the orthogonal 
or symplectic group and the associated dual symmetric space parametrized by antisymmetric or symmetric complex 
matrices, respectively [^, ^ 0. 

APPENDIX B: GENERALIZED COLOR-FLAVOR TRANSFORMATION 

Kaneko (see also Yan |Q) has derived the following remarkable integral identity: 

1 } ^ 

^-"W ^ 5.(A. + n. A.. y n {"y y.f' iy - ')") Mf)' 

1 r " 

= g„(y,,y,,l/O y n(^^'-^ani-^a)^ni-to.)^) {A^ixff' , (Bl) 
c 

where the constants are related as 

V, = ^^^±4^+N-2- V, = -^-N. (B2) 
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The normalization constant Skia, b, c) is given by the Selberg integral 



fe-i 



Sk{a,b,c) = Yi 



r(a + 1 + cj)Tib + 1 + cj)r(l + c{j + 1)) 
T{a + b + 2 + c{k + j - l))r(l + c) 



The integration contour C encircles the cut between Xa — and = 1. In our case (cf. Eq. (19)) 

/2 + n-l/2 



Ai.2 = ^1,: 



1 



h+h + l 
I 



N-2: 



V-2 



I 



N . 



(B3) 



(B4) 



We call the identity (Bl) fcrmionic replica "generalized color-flavor transformation". Indeed, for the RMEs values of 
the parameters (A = 1/2, 1, 2 and special values of A1.2) Eq. (Bl) essentially coincides with the fermionic replica version 
of Zirnbauer's color-flavor transformation pO| , p^ , after proper parametrization of the symmetric space elements and 
i nte gration out of the irrelevant angles. The questions whether there is a geometrical "dual pair" interpretation of Eq. 
(Bl) for arbitrary parameters, and whether there is a supersymmetric (say for rational A) or bosonic replica analogue, 
are currently open. 

In the large-iV limit we collect all the terms having the iV-th power into exp[— ^^^j^ S{xa,t)], where 



S{x, t) — — log(l — tx) — log 2: + log(l — x) . 



(B5) 



We then look for the stationary points of the "action" S{x) given by solutions of dxS = 0. A simple algebra gives for 
the stationary points 



x± = l±i 



1 - t 



t 



(B6) 



We then magnify the vicinity of t = by introducing 6 as t — sin^(0/ (2A^)) ~ 0"^ / (2NY and changing the integration 
variable x„ to (A„ as 



1 



1 - t 



(B7) 



The two saddle points (B6) are at (pa = 0,7r. Taking the limit N 00, one obtains for the action NS{xa,t) 
—i9 cos (j)a- A straightforward algebra yields 



(A„(a;) 



n 

a=l 



e'~ 



n 

.a>b 



sm 



1/A 



(B8) 



The closed contour C in the Xa-plane can be taken to be a circle, so that (f)a € [0,2tt]. As a result, one obtains Eq. 
( pT| ) of the main text. 
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